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For a brief time in its early evolution the Universe was a cosmic nuclear reactor. The expansion 
and cooling of the Universe limited this epoch to the first few minutes, allowing time for the 
synthesis in astrophysically interesting abundances of only the lightest nuclides (D, 3 He, 4 He, 
7 Li). For big bang nucleosynthesis (BBN) in the standard models of cosmology and particle 
physics (SBBN), the SBBN-predicted abundances depend on only one adjustable parameter, the 
baryon density parameter (the ratio by number of baryons (nucleons) to photons). The predicted 
and observed abundances of the relic light elements are reviewed, testing the internal consistency 
of primordial nucleosynthesis. The consistency of BBN is also explored by comparing the values 
of the cosmological parameters inferred from primordial nucleosynthesis for the standard model 
and for models with non-standard early Universe expansion rates with those determined from 
studies of the cosmic background radiation, which provides a snapshot of the Universe some 400 
thousand years after BBN ended. 
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1. Introduction 

During its early evolution the Universe passed through a very brief epoch during which its 
high temperature and particle density allowed for nuclear reactions among nucleons building the 
lightest nuclides (D, 3 He, He, 7 Li) in astrophysically interesting abundances. As a result the com- 
parison between the predicted and observed primordial abundances of the light nuclides provides 
a unique window on the early evolution of the Universe and a key probe of the standard models 
of cosmology and particle physics. Much later in the evolution of the Universe, after electrons 
combine with protons and alphas at "recombination", the cosmic background photons are free to 
propagate, forming the cosmic microwave background (CMB) radiation, the black body spectrum 
of radiation observed today at Tq = 2.725 K. Primordial, or big bang nucleosynthesis (BBN) and 
the CMB probe particle physics and cosmology at two very widely separated epochs in the evo- 
lution of the Universe: BBN when the Universe was only ~ 20 minutes old and the CMB some 
400 thousand years later. Comparisons of the predictions and observations of BBN and the CMB 
provide tests of the standard models of cosmology and particle physics and offer constraints on 
new physics. For a recent review of these tests and constraints and for further references, see [lit]. 

The standard model of cosmology uses General Relativity to predict the evolution of an 
expanding Universe filled with radiation, including three flavors of light neutrinos, and matter 
(baryons along with non-baryonic dark matter). Neither non-baryonic dark matter, dark energy (or 
a cosmological constant), or spatial curvature are relevant for our discussion here since they play 
no role in the physics and evolution of the early Universe. BBN in the standard model (SBBN) 
predicts the primordial abundances of the four light nuclides as a function of only one adjustable 
parameter, the "baryon abundance parameter", the ratio by number of baryons (nucleons) to CMB 
photons, T]b = n^/tiy = 10 10 T]io. T]b provides a measure of the baryon asymmetry of the Universe 
whose value is a goal of but, at present, not a prediction of models of particle physics. A large 
class of models of cosmology and/or particle physics which go beyond the standard models predict 
a non-standard expansion rate during the early, radiation-dominated (RD) Universe. The Hubble 
parameter, H, measures the universal expansion rate; for the early RD evolution the age of the 
Universe, t, and H are related by Ht = 1/2. A non-standard expansion rate may be parameterized 
by the "expansion rate parameter", 5 = H' /H, where 5 = 1 for the standard model. Historically, 
a non-standard particle content of the early Universe has been considered as the source of S ^ 1 
[01. In this case, S may be related to the "equivalent number of additional neutrinos", the "extra" 
energy density normalized to the energy density contributed by one flavor of standard model neu- 
trinos: AN V = (p' — p)/p v (N v = 3 + AN V ). S and AN V provide equivalent parameterizations of 
non-standard expansion rates. Prior to e annihilation, 5=1 + 7AN v /43 = 1 + 0.163AN V , while 
post-e 1 * 1 annihilation, S 2 = 1 +0.134AN V (see, e.g., However, keep in mind that AN V ^ 
need not be the result of extra neutrino flavors but could be due to other extensions of the stan- 
dard models. In particular, AN V need not be an integer and need not be positive. For example, if 
the gravitational constant were different in the early Universe (at BBN and/or at recombination) 
from its value today, G B bn/G = 1 + 0.163AN V and Grec/G = 1 + 0.134AN V . While SBBN only 
depends on rj^, BBN depends on S (or AN V ) as well. 

BBN at ~ 3 minutes and the CMB (recombination) some ~ 400 thousand years later provide 
complementary probes of the physics and early evolution of the Universe. This review of BBN ad- 
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dresses the following questions. Do the BBN-predicted abundances agree with the observationally- 
inferred primordial abundances (i.e., is BBN internally consistent)? Do the values of T]b inferred 
from BBN and the CMB agree? Do the values of AN V derived from BBN and the CMB agree and, 
are AAV(BBN) = AA^ V (CMB) = 0? 

2. SBBN-Predicted And Observationally-Inferred Primordial Abundances 

The observations of value in inferring the primordial abundances are diverse, from the oldest, 
most metal-poor stars in the Galaxy ( 7 Li), to Hn regions in the Galaxy ( 3 He) and extragalactic 
Hn regions ( 4 He), to cool, neutral (Hi) gas in the Lyman-alpha forest (D). To a greater or lesser 
extent, the material in these regions have experienced some post-BBN nuclear processing, poten- 
tially modifying the original relic abundances. The post-BBN evolution of deuterium is simple and 
monotonic: as gas is cycled through stars D is burned to 3 He and beyond As a result, the deu- 
terium abundance observed anywhere, at any time in the evolution of the Universe provides a lower 
bound to the primordial value and, the observed deuterium abundance should reveal a plateau at the 
primordial value in regions at high redshift and/or low metallicity where minimal stellar process- 
ing has occurred. In contrast, the post-BBN evolution of 3 He is considerably more complicated. 
While some of the 3 He incorporated into stars (along with any D which is burned to 3 He) survives 
in the cooler, outer regions, most is burned away in the hotter interiors and, newly synthesized 
3 He is produced by incomplete hydrogen burning in some stars. The net effect of these competi- 
tive processes depends on detailed models of stellar structure and evolution, along with models of 
galactic chemical evolution. If there is net production of 3 He in the course of chemical evolution, 
the 3 He abundance should correlate with the "metal" (i.e., CNO...) abundance which, in turn, is 
correlated with location in the Galaxy; for net destruction an anticorrelation would be expected. 
Since stars burn hydrogen to helium in the course of their evolution, the post-BBN evolution of 
4 He is simple and monotonic (similar to that of deuterium). The observationally-inferred abun- 
dance of 4 He should correlate with metallicity (i.e., CNO...), extrapolating to the primordial value 
at zero metallicity. Finally, like deuterium, most lithium is destroyed when gas is cycled through 
stars. However, some stars appear to be net producers of 7 Li and collisions between cosmic ray 
and interstellar nuclei also produce post-BBN lithium. All these effects must be kept in mind when 
using the observations to infer the primordial abundances. 

2.1 Deuterium 

Given its simple post-BBN evolution and the sensitivity of the BBN-predicted abundance to 
the baryon abundance parameter ((D/H)dp % L6 , so that a ~ 10% determination of (D/H)p yields 
a ~ 6% measurement of T]b), deuterium is the baryometer of choice. There is, unfortunately, a very 
small set of D abundances inferred from the spectra of only seven, high-redshift, low-metallicity, 
QSO Absorption Line Systems []|] (the observations require high resolution spectra on the world's 
largest telescopes). Even more problematic, the seven data points exhibit an unexpectedly large 
dispersion (given the quoted errors). These variations in D abundances show no obvious correlation 
with redshift or metallicity. From their data Pettini et al. [|5|] derive log = 5+log(D/H)p = 
0.45 ± 0.03, where the error has been inflated from the error in the mean to account for the large 
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Figure 1: The Izotov & Thuan 2010 (IT10) ||| helium and oxygen abundances. The solid line is the IT10 
best fit for a linear Y versus O/H relation; see the text. 

dispersion. In the context of SBBN this estimate of the primordial D abundance corresponds to a 
baryon abundance of 77i (SBBN) = 5.80 ±0.27. 

2.2 Helium-3 

For the baryon abundance parameter inferred from SBBN using the Pettini et al. deuterium 
abundance, the SBBN-predicted 3 He abundance is j 3P (SBBN) = 10 5 ( 3 He/H) P = 1 .08 ±0.04. Aside 
from the solar system [01, 3 He is observed in Galactic Hn regions via the spin-flip transition (the 
analog of the 21 cm line in neutral hydrogen) in singly-ionized 3 He. The most complete data set is 
that of Bania, Rood & Balser [J7|]. The gas in these Galactic H n regions has been processed through 
several generations of stars so, not surprisingly, the 3 He abundances reveal a large variation, likely 
due to post-BBN production/destruction. However, these 3 He abundances show no clear correlation 
with either metallicity {e.g., O/H) or location in the Galaxy (see, e.g., Figure 9 and the discussion 
in Steigman 2007 [jl|]). The lowest 3 He abundances, > 1.1 ±0.2, which are adopted by Bania, 
Rood & Balser as an estimate of (or, a lower bound to) the primordial abundance, are in excellent 
agreement with the SBBN-predicted abundance, providing support for the internal consistency of 
SBBN (two predicted abundances, one free parameter). However, the lack of a correlation between 
the 3 He abundances and metallicity or location in the Galaxy is puzzling. 

2.3 Helium-4 

As the second most abundant nuclide in the Universe 4 He is observed in many astrophysical 
environments. The most useful observations for BBN are those of the hydrogen and helium recom- 
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Figure 2: The lithium abundances ([Li] = 12+log(Li/H)) as a function of metallicity ([Fe/H]) as derived 
from observations of very metal-poor stars in the Galaxy. See the text for references. The horizontal band 
shows the ± 1 a SBBN prediction. 

bination lines from low-metallicity, extragalactic Hn regions. This is a field with a long history; 
see [jl|] for references. With the presently available large data set it has become increasingly clear 
that systematic uncertainties (which are always difficult to quantify) dominate the statistical errors. 
Recently, Izotov & Thuan [Jsj] used their data of 93 spectra for 86 low-metallicity extragalactic 
H ii regions to derive the helium and oxygen abundances shown in Figure [j] A linear fit to the Y 
- O/H data extrapolated to zero metallicity yields Y P = 0.2565 ± 0.0010(s/af) ± 0.0050(jyjf) [§]. 
For the comparison here with BBN I combine the errors linearly, adopting Yp = 0.2565 ±0.0060. 
These results are consistent with the recent determinations of Yp from a much more limited data set 
(< 9 H ii regions) by Aver, Olive & Skillman [Q]. For SBBN with the D-determined baryon density 
parameter, Yp(SBBN) = 0.2482 ± 0.0007. While the central value of the SBBN-predicted abun- 
dance appears to be in conflict with the observationally-inferred value, given the relatively large 
uncertainty in the latter estimate, the difference is only at the ~ 1.4a level. SBBN is consistent 
with the observationally-inferred primordial abundances of D, 3 He, and 4 He (within the errors). 

2.4 Lithium-7 

For the D-determined baryon density parameter the SBBN-predicted lithium abundance is 
[Li] sbbn = 12+log(Li/H)sBBN = 2.65 ±0.06. In Figure ^| are shown the lithium abundances de- 
rived from observations of the oldest, most metal-poor stars in the Galaxy [ |T0| , 11, 12, 13]. These 



stars should provide a sample of the "lithium plateau", the nearly primordial abundance in systems 
which have experienced very little stellar processing. It is difficult to identify a lithium plateau 
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from the data shown in Figure |2| and it is clear that none of the observationally-inferrred lithium 
abundances is even close to the SBBN-predicted value. The difference is ~ 0.5 — 0.6 dex or, a 
factor of 3 - 4. Lithium is a problem. However, since the target stars have had > 10 Gyr to modify 
their surface abundances by mixing with the lithium-depleted interior material, it is unclear if the 
problem is one for cosmology or particle physics or, for stellar astrophysics. 

2.5 Summary for SBBN 

The observationally-inferred primordial abundances of the light nuclides are subject to diffi- 
cult to quantify systematic uncertainties along with the usual statistical errors. The fact that D, 
3 He, 4 He, and 7 Li are observed by different astronomical techniques in different targets ensures 
that their derived abundances are not affected by the same systematics. Within the errors, the 
SBBN-predicted and observationally-inferred primordial abundances of D, 3 He, and 4 He are in 
good agreement, providing support for the standard models of particle physics and cosmology 
(e.g., AAV = 0). However, lithium is a problem. Since non-standard BBN offers a second free 
parameter, S or AAV, it is of interest to see how the BBN abundances constrain the combination of 
tjb and AAV and, if this additional freedom can help to ameliorate the lithium problem. 



3. Non-Standard BBN 



For the more general BBN case the relic abundances depend on two adjustable parameters, the 
baryon density parameter T7b and the expansion rate parameter S (or, equivalently, AAV). m Figure 
|3| are shown the D and 4 He isoabundance contours, vdp (dotted, blue) and Yp (solid, red) in the 
S — TJio plane Jl4| , |l|]. The point with the error bars represents the observationally-inferred D and 
4 He abundances from §||. Since the vdp and Yp contours form a grid in the S — T]io plane, the D 
and 4 He abundances constrain S (or, AAV) and Tjio, as may be seen from Figure [|. The adopted D 
and 4 He abundances predict Tjio(BBN) = 6.07 ±0.33 and AAV = 0.62 ± 0.46 (= at ~ 1.3a). 

Also shown in Figure [3] (solid, green) are the 7 Li isoabundance contours. The D and 4 He 
constrained, BBN-predicted lithium abundance with both T7b and S (AAV) as free parameters yields 
[Li]BBN = 2.66 ± 0.06. Lithium is still a problem. 



4. Comparison Between BBN And The CMB 

The CMB offers a probe of the later, but still early, evolution of the Universe which comple- 
ments that from BBN. In particular, the CMB temperature anisotropy spectrum is sensitive to the 
values at recombination of both the baryon abundance parameter and the expansion rate parameter 
(or, AAV). 



From the 7-year WMAP data Komatsu et al []15|] find r\ m = 6. 190 ±0.145. If it is assumed that 



AV = 3 (SBBN), for which T7i (SBBN) = 5.80 ±0.27, then SBBN and the CMB agree to within 
~ 1.3a on the value of the baryon density parameter at a few minutes and some 400 thousand years 
later. The comparison of the SBBN and CMB inferred values of T]io are shown in the likelihood 
distributions of Figure ||. However, Komatsu et al find some evidence (at the ~ 1.5a level) in 
support of AAV + 0. For BBN (AAV + 0), T]io(BBN) = 6.07 ±0.33, which is in even better 
agreement with CMB result for Tjio, as may be seen in Figure |5[ 
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Figure 3: The BBN isoabundance contours in the S versus 7710 plane |[l4| [jj for y DP = 10 5 (D/H) (dashed, 
blue) and Y (solid, red), vdp decreases from left to right from 4.0 to 3.0 to 2.0 respectively. Y increases 
from bottom to top, from 0.24 to 0.25 to 0.26 respectively. Also shown are the isoabundance contours (solid, 
green) for [Li] = 12+log(Li/H). [Li] increases from left to right from 2.6 to 2.7 to 2.8 respectively. The data 
point with error bars corresponds to the adopted D and 4 He abundances; see the text. 



For BBN and the observationally-inferred D and 4 He relic abundances, AN V = 0.62 ± 0.46, 
which is consistent with AN V = at ~ 1.3a. From the imprint on the CMB at recombination 
Komatsu et al [jl|] find AN V = 1.30 ±0.87 (see also [||]), which differs from AN V = by ~ 1.5a. 
The very good agreement between BBN and the CMB, A/Y V (CMB) - A/V V (BBN) = 0.68 ± 0.98, as 
well as their overlap with AN V = 0, is illustrated by the likelihood distributions in Figure 0. 



5. Discussion 

The predictions of SBBN (AN V = 0) are consistent, within the errors, with the observationally- 
inferred relic abundances of D, 3 He, and 4 He, as well as with the value of the baryon density 
parameter derived from the CMB. But, lithium is a problem since its SBBN predicted abundance 
differs from that derived from observations of very metal-poor halo and globular cluster stars by 
a factor of three or more. The problem could be in the stars (e.g., depletion or dilution of surface 
lithium during the lifetime of the oldest stars in the Galaxy) or in the cosmology (e.g., late-time or 
renewed BBN initiated by the decay of long-lived massive particles). 

For non-standard BBN with AN V 7^ comparison between the predicted and observed abun- 
dances of D and 4 He constrains the allowed values of the baryon density parameter (t]b) and the 
expansion rate parameter (S) when the Universe was only a few minutes old. The BBN results 
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Figure 6: The likelihood distributions for N v = 3 + AN V inferred from BBN (red) and the CMB (blue). The 
vertical line (magenta) is for the standard model at N v = 3. 

are consistent with the CMB inferred values of these parameters at recombination, some 400 thou- 
sand years later. But, even in this case lithium remains a problem. Lithium aside, the very good 
agreement between these independent probes of cosmology during widely separated epochs in the 
evolution of the Universe validates using the combined constraints from BBN and the CMB to 
probe non-standard models of cosmology and particle physics. 

5.1 Comparing The Universe At BBN And At Recombination 

During most epochs in the evolution of the Universe entropy is conserved. The number of 
CMB photons in a comoving volume provides a measure of the entropy in that comoving volume. 
A comparison of the number of photons in the comoving volume at BBN and at recombination 
provides a test of entropy conservation. Of course, it is necessary to define the size of the comoving 
volume, which can be accomplished by identifying the number of baryons (nucleons) in it, so that 
N r = Nb/tjb- Given the very good constraints on baryon non-conservation, the number of baryons 
in a comoving volume is (should be!) preserved from BBN to recombination, so that a comparison 
of tjb derived from BBN and from the CMB constrains any entropy production in the intervening 
epochs: N 7 (REC)/N 7 (SBBN) = 0.94 ± 0.05 and N r (REC)/N r (BBN) = 0.98 ± 0.06. 

Similarly, the combined results of BBN and the CMB may be used to compare the value of the 
gravitational constant, G, at BBN and at recombination with each other and with the present value 
(G ). Since G B bn/G = 1 +0.163AN V (BBN) and G RE c/G = 1 +0.134AN V (REC) [§, the values 
of AAV at BBN and at recombination constrain the early Universe strength of gravity, leading to: 
Grec/Gbbn = 1.07 ±0.13, Gbbn/Gq = 1.10 ±0.07 and, Grec/G = 1.17 ±0.12. 



9 



Primordial Nucleosynthesis 



Gary Steigman 



As a last example of the value of combining the constraints from BBN and the CMB con- 
sider the effect of a massive particle which decays after BBN but before recombination. If the 
decay occurs too late for the relativistic decay products to be thermalized (producing extra CMB 
photons whose presence has already been constrained above), the energy they carry will, never- 
theless, contribute to the total energy density when the Universe is radiation dominated. Since 
(Pr/Pr)bbn = 1 +0.163AN V (BBN) and (p R /p R ) RE c = 1 +0.134AN V (REC), AN V measures the 
presence of "extra" relativistic energy. Comparing AN V at BBN with its value at recombination it 
is seen that they agree, (p r /Pr)bbn = (p r /Pr)rec within 0.5a. 

6. Conclusions 

For N v R3 3 BBN is internally consistent (but lithium is a problem!) and is in agreement with 
the CMB. The very good agreement between BBN and the CMB permits their use in combination 
to constrain some examples of non-standard models of cosmology and particle physics. While 
celebrating the success of BBN it should be kept in mind that some challenges remain. For example, 
why is the dispersion among the observed deuterium abundances so large? Or, why don't the 3 He 
abundances observed in Galactic H n regions correlate with the oxygen abundances or with location 
in the Galaxy? And, for 4 He, the second most abundant element in the Universe, how large are 
the systematic errors in its observationally-inferred primordial abundance and, are there observing 
strategies to reduce or eliminate at least some of them? This active research area of importance for 
cosmology and particle physics would benefit greatly from more data. 
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